Homopolar oscillating-disc dynamo driven by parametric resonance 
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We use a simple model of BuUard-type disc dynamo, in which the disc rotation rate is subject 
to harmonic oscillations, to analyze the generation of magnetic field by the parametric resonance 
mechanism. The problem is governed by a damped Mathieu equation. The Floquet exponents, which 
define the magnetic field growth rates, are calculated depending on the amplitude and frequency of 
the oscillations. Firstly, we show that the dynamo can be excited at significantly subcritical disc 
rotation rates when the latter is subject to harmonic oscillations with a certain frequency. Secondly, 
at supercritical mean rotation rates, the dynamo can also be suppressed but only in narrow frequency 
bands and at sufficiently large oscillation amplitudes. 
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In dynamo experiments, a high driving power is neces- 
sary to achieve the self-excitation of the magnetic field. 
The ensuing liquid metal fiow is usually strongly turbu- 
lent. In both Riga and Karlsruhe dynamo experiments 
the turbulent fluctuations were partly inhibited by 
the internal walls, whereas in the Cadarache experiment 
01 , the absence of such walls resulted in large-scale flow 
fluctuations The effect of flow fluctuations on the 
dynamo threshold has been addressed in several recent 
studies 0, i, 0, S, i, M, [HI, E, H [il • Solving the kine- 
matic dynamo problem for a given non-stationary flow 
usually governed by the Navier-Stokes equations shows 
that turbulence generally has an adverse effect on the 
dynamo excitation, unless the fluctuations are strong 
enough to drive the dynamo by themselves without any 
mean flow. In the latter case, we speak of a fluctua- 
tion dynamo [isl . whose experimental implementa- 
tion seems hardly feasible because of the high excitation 
threshold. However the possibility that fluctuations ex- 
cite the magnetic fleld by the parametric resonance mech- 
anism can not be excluded. Parametric resonance 
has been proposed in the somewhat different context of 
spiral galaxies as promoter of bisymmetric magnetic fleld 
structure [H, HH . 

In this paper, we use a simple model of the Bullard- 
type disc dynamo |23| to show that the magnetic fleld 
can indeed be excited by the parametric resonance mech- 
anism, even when relatively small harmonic oscillations 
are added to signiflcantly subcritical disc rotation rates. 

Consider a Bullard-type disc dynamo [131 which con- 
sists of a solid conducting disc rotating with a generally 
time-dependent angular velocity z{t) about its axis, and 
a wire twisted around the axle and connected by sliding 
contacts to the rim of the disc and the axle as shown in 
Fig. [TJ The disc is assumed to be segmented so that 
azimuthal current can flow only at its rim. This corre- 
sponds to the modiflcation of the Bullard disc dynamo 
suggested by Moffatt in order to eliminate exponential 
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Figure 1: Sketch of a homopolar disc dynamo. 



growth of the magnetic fleld in the limit of a perfectly 
conducting disc [23|- The system is described by the fol- 
lowing set of dimensionless equations (for details see Ref. 



X = r{y-x), 

y — xz + mx — {m + l)y, 

z — g[l + x{mx — (m -f- l)y)] 



kz, 



(1) 



where x and y are magnetic fluxes through the loop made 
by the wire and the rim of disc, respectively; z is the 
dimensionless angular velocity of the disc; r accounts for 
the resistance of the disc relative to that of the loop, and 
m characterizes the relative mutual inductance of the disc 
and the loop; the dot stands for the time-derivative d/dt. 
The disc is driven by a generally time-dependent torque 
g, and braked by a viscous-type friction characterized 
by the coefficient k which is necessary for the structural 
stability of the system [131 ■ Henceforth, we assume the 
friction to be strong with respect to the inertia of the 
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Figure 2: Growth rate 7 versus the reciprocal frequency (a) 
marginal mean rotation rate at ao = 0. 

disc accounted for by i in Eq. which, thus, results 
in z = zq [1 + x{mx — (m + l)y)] , where zq = g/k. The 
remaining two Ist-order ODEs in ^ can be combined 
into a single 2nd-order DufRng-type equation [3] with a 
non-linear friction 

x+{l + I3x^)x-ax + \x^ = Q, (2) 

where x and t are rescaled by (m + 1 + r) and (m + 
1 + r)^"^ , respectively, and a = r(zo — l)/(m + 1 + r)^, 
/? = zo(to + 1)(to + 1 + r), and A = rzo- Further, we 
focus on the evolution of small initial perturbations of 
the magnetic field characterized by x ^ 1, for which Eq. 
([2]) can linearized by setting (3 = \ — Q. Then the only 
remaining parameter a depends directly on the deviation 
of the disc rotation rate from its critical value a = 0. 
For a > 0, a small initial magnetic field starts to grow 
exponentially provided that the disc rotates steadily |22l |. 
In this study, we are interested in how the generation of 
the magnetic field is affected by the unsteadiness of the 
disc rotation 

a — + 6 cos{ijjt) , (3) 

which besides the mean part ao contains also an oscil- 
latory component with the amplitude 5 and the circular 
frequency u. Then the linearized Eq. ^ reduces to a 
damped Mathieu equation 

X + X ~ (ofQ + 6 cos{ujt))x ~ 0. 

Using the substitution x{t) = exp{—t/2)x{cdt/2), the 
equation above can be transformed into the canonical 
Mathieu equation 

X+[a-2gcos(2T)]x = 0, (4) 

where a — —(1 + 4Q!o)/a;^, q = 26/lo'^ and r ~ ut/2. 
According to Floquet theory, a particular solution to Eq. 
(0]) can be written as x{t) = exp(wr)/(T), where /(r) 
is a TT-periodic function and v is the Floquet exponent — 
both dependent on the parameters a and q. According 
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and the critical amplitude 5c versus the frequency (b) for the 



to this solution, the amplitude of the magnetic field x{t) 
evolves exponentially in time with the maximum growth 
rate 7 = — 1) /2, where the modulus accounts 

for the time-refiection symmetry of Eq. (0]) [2^. Thus, 
the amplitude of the magnetic field grows exponentially 
when 7 > 0, whilst the marginal state is defined by 7 = 0. 
We use the Maple computer algebra software to calculate 
Floquet exponent which defines the growth rate 7 de- 
pending on the amplitude 5 and the frequency u. Next, 
we find the critical oscillation amplitude 5c as the func- 
tion of frequency u) for fixed values of the mean rotation 
Qfo by solving numerically equation = 1 corre- 

sponding to 7 = 0. Eventually, we determine the mini- 
mal oscillation amplitude Jmin and the corresponding fre- 
quency at which an exponentially growing magnetic field 
first appears. The corresponding numerical results are 
presented and discussed below. 

We start with a marginal mean disc rotation rate 
ao = 0, which corresponds to the dynamo excitation 
threshold when the disc rotates steadily, i.e., 7 = when 
(5 = 0. As seen in Fig. [2]Ja), the disc oscillations about 
the critical rotation rate with a sufficiently small ampli- 
tude {5 = 0.25) brings the growth rate to a constant 
level below zero as the frequency is reduced (recipro- 
cal frequency increased). As the oscillation amplitude 
increases, first, the growth rate splits into separate fre- 
quency bands whose width decreases as ~ I/lo for 0. 
In order to show this increasingly fine-scale structure of 
the growth rate as — > 0, we use the reciprocal frequency 
which is proportional to the period of disc oscillations. 
Secondly, the growth starts to increase with the oscilla- 
tion amplitude and approaches zero again at (5 w 0.75 for 
a certain critical frequency. Further increase in the oscil- 
lation amplitude to (5 = 1 results in the appearance of sev- 
eral frequency bands with positive growth rates (7 > 0) 
which are shown as filled regions in Fig. Wis). The criti- 
cal amplitude 5c at which the growth rate turns zero, is 
shown for ao = in[2]Jb) against the oscillation frequency. 
Marginal state with 7 = 6 corresponds to the boundaries 
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Figure 3: Growth rate 7 versus the reciprocal frequency (a) an 
subcritical mean rotation rate at qq = —10. 

of the separate frequency bands shown by different col- 
ors. Growth rate is positive corresponding to the dynamo 
action inside the filled frequency bands which approach 
each other closely as the oscillation amplitude increases. 
Note that only a certain number of first instability bands 
are shown in Fig. ^h). There is an infinite sequence of 
similar instability bands of decreasing width as w ^ 0. 
Thus, the range of unstable frequencies, which is inter- 
vened by infinitely many, much narrower stability bands, 
extends down to w = 0, whereas it is bounded from above 
by the first instability band. The critical oscillation am- 
plitude, which is the lowest for the first instability band, 
rises to an asymptotic value depending on a as w ^ 0. 
The minimal value of the critical oscillation amplitude 
and the corresponding frequency at which it occurs are 
plotted in Fig. [4] versus subcritical (negative) values of 
a. 

The reduction of disc rotation rate to a moderately 
subcritical value of ao = —1, results in the increase of 
the minimal oscillation amplitude necessary for the dy- 
namo action (7 > 0) up to (5c ~ 2. It means that the 
maximum disc rotation rate ^ during the oscillation 
cycle temporally exceeds the critical value a = for a 
steady rotation. This, however, changes when the disc 
rotation rate is reduced further down to ao = —10, for 
which the growth rates and the corresponding genera- 
tion bands are shown in Fig. [3l In this case, a positive 
growth rate band is seen in Fig. El^a) already at Sc = 8. 
This implies that a < during the whole cycle of disc 
oscillation. Thus, the dynamo appears at a maximum 
disc rotation rate below its critical value for a steady ro- 
tation. There is a range of frequencies seen Fig. Eljb), at 
which a subcritical growth of the magnetic field is possi- 
ble with the the oscillation amplitude Sc < —ao = 10. As 
seen in Fig. lU such a subcritical excitation of magnetic 
field appears already at ao ^ —4, where the minimal re- 
quired oscillation amplitude (Smin becomes smaller than 
—ac, which is shown by the dashed fine. Note that for 
a strongly subcritical ao and steadily rotating disc, an 
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initial perturbation of the magnetic field oscillates with 
the period 0(1/^/— ao) which is much shorter than the 
characteristic damping time 0(1). Thus, for such ao the 
damping of magnetic field oscillations becomes relatively 
slow that enables parametric excitation of the magnetic 
field by a relatively weak modulation of the disc rotation 
rate. 

For ao < —0.25, the bands of unstable frequencies are 
associated with the subharmonics of Mathieu equation 
defined by a^/^ = 1,2,3, ... in Eq. ^ Thus, the 

most unstable band corresponds to the first subharmonic 
with LUi = y^—iao — 1 which is seen in Fig. |4]to approx- 
imate oJc well for u! > 1. In addition, it is interesting to 
note that the critical frequency, ujc, at which i5min occurs, 
is numerically very close to (Smin itself in the whole range 
of ao < 0. 

For a supercritical disc rotation rate with ao = 1, the 
growth rate 7, which is seen in Fig. [DJa) to be positive 
for a steadily rotating disc {S — 0), first reduces as the 
oscillation amplitude is increased to S — 2. Moreover, 7 




Figure 4: Minimal oscillation amplitude and critical frequency 
versus subcritical rotation rate parameter (a < 0). 
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is seen to become negative within certain, relatively nar- 
row frequency bands. These negative growth rate bands 
(7 < 0), where the dynamo is suppressed by the disc 
oscillations, are shown in Fig. ^a) by filled curves. As 
seen in Fig. IMp), the width of the first suppressed fre- 
quency band noticeably increases while the whole band 
shifts towards the high-frequency range, i.e., low recip- 
rocal frequencies, as the oscillation amplitude increases. 
At the same time, the width of the subsequent frequency 
bands, where the dynamo is suppressed, remains very 



small, while the spacing between the adjacent suppres- 
sion bands reduces in terms of the reciprocal frequency. 

In conclusion, we have shown in this study that dy- 
namo can be excited by the parametric resonance at 
considerably low velocities of the conducting medium 
when the latter is subject to harmonic oscillations. 
Even though this was demonstrated for a simple model 
Bullard-type disc dynamo, we expect that similar mech- 
anism may also be relevant for more realistic dynamos. 
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